Recoil-induced lasing 
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The interaction of an atomic gas confined inside a cavity with a strong electromagnetic field is 
numerically and theoretically investigated in a regime where recoil effects are not negligible. The 
spontaneous appearance of a density grating (atomic bunching) accompanied by the onset of coherent 
backpropagation is found to be ruled by a continuous phase-transition. Numerical tests allow us 
to convincingly prove that the transition is steered by the appearence of a periodic atomic density 
modulation. Consideration of different experimental relaxation mechanisms induces us to analyze 
the problem in nearly analytic form, in the large detuning limit, using both a Vlasov approach and 
a Fokker-Planck description. The application of our predictions to recent experimental findings, 
reported in [Phys. Rev. Lett. 91, 18360f (2003)], yields a semi-quantitative agreement with the 
observations. 
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I. INTRODUCTION 

Recoil in the interaction between atoms and electro- 
magnetic fields is almost exclusively associated with the 
idea that momentum transfer is a way of slowing down 
atoms to extremely low-temperatures However, non 
trivial consequences on light propagation have been un- 
covered as well. For instance, in the realm of single-atom 
effects, the so-called Recoil Induced Resonances (RIR) 
have been investigated both theoretically and exper- 
imentally |3j. Moreover, it has been conjectured that 
atomic recoil can collectively give rise to coherent propa- 
gation through backreflection from a spontaneously gen- 
erated density grating: the so called Collective Atomic 
Recoil Laser (CARL) Finally, connections between 
the two approaches have been investigated in p| . 

The very first observation of spontaneous backward 
emission was made in a ring cavity experiment 6] , where 
the very low "transverse temperature" might have been 
responsible for CARL amplification. More detailed ex- 
perimental accounts have been reported in Refs. @,|^], 
although it was not possible to establish whether the 
generation of the backward field was due to the spon- 
taneous formation of a density grating, or viceversa. In 
fact, although it was later ascertained that recoil plays 
a prominent role [f| in the underlying physics Q, ; an 
alternative interpretation based on the formation of a 
standard polarization grating could not be entirely ruled 
out |lfj| . Only recently has the first convincing evidence 
of CARL been given in a beautiful experiment performed 
on a sample of cold Rubidium atoms \H\ . 

On the theoretical side, too, no final conclusion could 
be drawn on the CARL features, because the first model 
allowed only investigating the transient regime. Station- 
ary states have been obtained for the first time by includ- 



ing collisions with an external buffer gas By simu- 
lating a low-temperature (a few mK) sodium vapour in 
the detuning range previously considered in the (numeri- 
cal) literature, a noncquilibrium phase transition leading 
to a stationary nonzero backward field (above a given 
threshold) was therein identified. However, such a col- 
lective phenomenon could not be linked to the onset of 
a density grating, but rather to the birth of a nontrivial 
polarization grating [T^ . 

The overall scenario has been observed also in a more 
general model accounting for the input-field dynam- 
ics [l3| and has been successively clarified in [l4( with 
the help of the elimination of the atomic variables; in- 
deed, this step rendered the derivation of an effective 
free-energy potential for the backward field amplitude 
possible, thereby clearly showing the existence of a phase 
transition. 

In this paper we show numerically and theoretically 
that this latter model also accounts for the CARL and al- 
lows one to understand its onset in terms of a continuous 
phase transition characterized by two effective variables. 
The resulting estimate of the critical point is compati- 
ble with the experimental results reported in The 
results presented here provide the first substantial and 
systematic progress towards a proof of the existence, in 
stable form, of a collective interaction, which, predicted 
years ago |J, had so far proven to be elusive. 

More precisely, we investigate the transition which 
gives rise, spontaneously - i.e., without an external probe 
field -, to a coherent electromagnetic field counterpropa- 
gating with respect to the injected one. We choose to con- 
sider the ring-cavity-based model introduced in [l3L ITo| 
for two reasons. On the one hand, the available single- 
pass model has intrinsic shortcomings 0] which may 
be of little consequence in particular circumstances, but 
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which renders it nonetheless less attractive. On the other 
hand, the most recent experimental results [Tlj . which 
can be considered at the present time as the first and 
only available observation of a true CARL system, have 
been obtained in a ring cavity. Hence, the choice of a 
bidirectional ring cavity model imposes itself. 

A strong concern, which tainted all previous predic- 
tions and experimental attempts at studying this prob- 
lem, regards the influence of external perturbations (col- 
lisions, finite transit time in the interaction volume, "op- 
tical friction", etc.) on the collective dynamics. One can 
conjecture that their presence could destroy the phase 
transition or, at least, affect it so strongly as to shift 
it to an unphysical region of parameter values. Our pur- 
pose in this work is to investigate such mechanisms using 
different approaches and mathematical/numerical tools, 
to obtain a convincing picture. Section [n] introduces 
the model we used in the investigation and discusses, in 
subsectio rlll Al the main physical relaxation mechanisms 
that may take place in an experiment. The appropriate 
modeling tools for analyzing the two main mechanisms 
are introduced here. Section ITTT1 is devoted to a system- 
atic numerical analysis which allows us to highlight the 
origin of the forces that give rise to a transition of the 
CARL-type - i.e., caused by the appearence of a spatial 
density grating - as opposed to the one recently pre- 
dicted 0, 0, based on an unusual polarization grat- 
ing and possibly responsible for high-temperature ex- 
perimental observations 0, 13 • The numerical approach 
makes it possible to isolate the different source terms 
and identify in which range of temperature each of them 
is predominant. Analytical descriptions present the ad- 
vantage of offering a better understanding of the physics 
driving the system. Since the full problem is obviously 
analytically unmanageable, it is useful to study some lim- 
iting situations. This is done in Section lTVI where we con- 
sider the case of large detuning between fields and atomic 
resonance. In this spirit, we describe the perturbations 
introduced by transit broadening (subsection IIV A|) and 
show that a second-order transition appears which can 
be described by expanding the spatial inhomogeneities in 
Fourier modes. The first and most relevant mode coin- 
cides with the bunching parameter, previously introduced 
to quantify the density grating |4| . Losses introduced by 
an optical molasses are studied in subsection IIV Bl with 
the help of a Fokker-Planck equation that can be ana- 
lytically solved in the large dissipation limit. Predictions 
resulting from the preceding analyses are offered in Sec- 
tion[3 together with a comparison with the recent exper- 
imental results [TTf. Some conclusions and perspectives 
are presented in Section IVT1 



II. MODEL 

The starting model involves four dimensionless vari- 
ables (one being complex) accounting for the single atom 
dynamics plus two complex equations describing the 



amplitude of forward (xf) and backward (x b ) uniform 
fields [L| 

e 3 = pj, 

Pj = -$te[(x f e ie * -x b e- ie *)s*](\A a \/2G), 

Sj = G{x f e i9 i +x b e-^)d 3 -(T ± -iA a ) Sl , (1) 

dj = -G$te[(x f e ie * +x b e- i6 i)s*]-T ]l (d j + l), 

x f = -{l + iA c )x f +Y + C(sc- ie ), 

x b = -(l + iA c )x b + C(se l9 ), 

where 1 < j < N and the angular brackets denote an av- 
erage over the atomic ensemble. Time and related param- 
eters are expressed in units of the photon life-time inside 
the cavity, namely: the decay rates of the atomic po- 
larization (rj_) and population inversion (Tu), the detun- 
ing A a between input field and atomic frequencies, and 
the cavity detuning A c (i.e. the distance from the near- 
est cavity resonance). Moreover, G = ^/mn\ A a \/h/k 
and C = aT ±L / (GT) are the two coupling constants (m 
is the atomic mass, k the field wavenumber, a the unsat- 
urated absorption rate per unit length, L the length of 
the atomic medium, and T the cavity's transmittivity). 
The amplitude Y of the injected field is scaled such that 
Y 2 = AV 2 V/{AG 2 T K 2 h 2 ), where V is the atomic dipole 
moment, V is the input power, and A the beam sur- 
face, expressed in physical units. An equivalent scaling 
is used for the forward and backward intracavity com- 
ponents, Xf and x b . Finally: (i) 9j = kzj is the scaled 
atomic position; (ii) the atomic momentum Pj is rescaled 
to pj = kPj I '(run); (iii) Sj and dj are the atomic polari- 
sation and population inversion, the first being expressed 
in units of T>. 



A. Relaxation mechanisms 

While dissipation mechanisms affecting the internal de- 
grees of freedom and the field amplitudes are accounted 
for by the above model, no losses acting on the atomic 
momentum are therein included. As a result, atoms 
are continuously accelerated by radiation pressure, and 
Eqs. can thereby describe only the transient regime 
corresponding to the early stages of CARL amplification. 

In actual experimental situations, relaxation processes 
may be induced by various mechanisms, the main ones 
being: (a) collisions with a buffer gas (those between 
optically active atoms can be neglected under all realis- 
tic circumstances), (b) finite residence time in the in- 
teraction volume (added to some relaxation occurring 
outside it), and (c) "viscous" losses (e.g., as in an op- 
tical molasses) . Collisions with a buffer gas require fairly 
high pressures and have been often used in the inves- 
tigation of nonlinear coupling between atoms and light 
(cf., e.g., [TtI [l8|L In a numerical scheme, this relax- 
ation processes can be implemented as a coupling with a 
thermal bath, efficiently simulated with "molecular dy- 
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namics" [l9j techniques. More precisely, a random se- 
quence of inter-collision times is independently gener- 
ated for each atom according to an exponential distribu- 
tion with average value equal to 1/7 C . At each collision, 
the momentum of the colliding, j-th, atom is randomly 
reset according to the Gaussian distribution Qo{pj) — 
exp {-p 2 /(2T)} /VZttT where T = k B fk 2 /(mn 2 ) is the 
rescaled temperature of the buffer gas (fee is Boltzmann's 
constant, T is the temperature in physical units); more- 
over, the phase of the atomic polarization Sj is also re- 
set to a value uniformly distributed in the whole range 
[0, 2ir], as \sj\ remains unchanged. Notice that this algo- 
rithm amounts to implicitely assuming that active and 
buffer-gas atoms have equal masses. 

The finite transit time (b) introduces effective losses 
and has certainly played the main role in several exper- 
iments 0,00- In sucn a situation, the atoms spend a 
more or less small fraction of time inside the interaction 
volume, which they enter and exit at random positions 
and in a random state (both for the velocity and for the 
internal degrees of freedom) . Outside the interaction re- 
gion, the atoms suffer thermalization, either through col- 
lisions with other atoms (e.g., a low pressure buffer gas), 
with a container, or are simply renewed through "fresh" 
atoms coming from a beam; hence, the effective relax- 
ation comes from the loss of "synchronized" atoms (both 
for the external and the internal degrees of freedom) com- 
pensated, in average, by the entrance into the interaction 
volume, of "random" atoms. This phenomenon can be 
reinterpreted as a "reset" in the external and internal 
degrees of freedom of the atom in question. At the mi- 
croscopic level, one can still adopt a description based on 
the thermalization scheme discussed in (a) , by modelling 
each random exit/entrance in the active volume as a "col- 
lision" ; the main difference with the previous case is that 
all atomic variables have now to be reset. At the macro- 
scopic level, the mathematical tool most suited for such 
a description is a Vlasov equation, where the probability 
distribution for the atomic position and momentum are 
dynamically investigated. These considerations are the 
basis for the discussion of Section IfV Al 

The third mechanism (c), recently employed in 
is an optical molasses, added in the experiment to ob- 
tain steady state backward emission. Here, the "molasses 
field" slows down the atoms by the standard cooling ac- 
tion while introducing a "friction" which prevents the 
Doppler shift from growing so large as to render the in- 
teraction with the field entirely negligible, ft has been 
shown [20| that under the assumption of a small satu- 
ration parameter of the molasses field, the momentum 
dynamics is well described by a Langevin equation (com- 
pare with the second of Eqs. Q), 

pj = $te[(x b e- ie * -x f e ie *)s*] (\A a \/2G) -<ypj +r ]j , (2) 

where rjj is a standard white noise, 

( Vj (t)) = 0, (3) 
(Vj(h)Vk(t 2 )) = 2jTS jk 6{t 1 -t 2 ), 



(5jk is the Kronecker symbol) accounting for the in- 
evitable fluctuations that thcrmalize the atom to a finite 
temperature T. In Sect ion Ft V Bl we will study this model 
in the limit of large detuning, a condition well verified in 
the experiment 



III. NUMERICAL RESULTS 

We now intend to explore the predictions of the model, 
Eqs. n from a purely numerical point of view to pro- 
vide a starting point for the later analysis. In particular, 
we are interested in studying the existence of a steady 
state bunching in the presence of a backward propagat- 
ing field. We remark that the observation of a stable 
value of bunching in a bidirectional ring cavity has never 
been reported. Previous observations, obtained without 
cavity 4], have shown a strongly oscillating, irregular be- 
haviour for this variable. In the other kind of observed 
transition \l2\ , instead, the 8 distribution remains 
flat (within the fluctuations imposed by the sample size) 
both below and above threshold. Hence, the mere numer- 
ical observation of a stable, steady bunching represents 
an important step towards generalizing the CARL pre- 
dictions [j]. 

Such an observation is, however, not sufficient in it- 
self. The question arises naturally, whether the backward 
field is originated by the growth of a nonhomogcnous 
longitudinal atomic distribution, or whether the spatial 
inhomogeneity results as a consequence of a preexisting 
optical standing wave, generated by other source terms. 
Although apparently obvious, this point is far from being 
trivial. Indeed, any standing wave is going to generate a 
whole hierarchy of gratings (atomic polarization and in- 
version, index of refraction and - if the atoms are free to 
move - atomic density) , but one needs to identify the true 
source. In this system, a source term which gives rise to 
a collective state has already been identified in the non- 
trivial grating in the atomic polarization's phase (l2lll4| . 
while a standard polarization modulation (without phase 
transition) has been proposed as a possible mecha- 
nism for interpreting experimental observations 0, Q • 

In order to characterize the appearence of a spatial 
inhomogeneity in the atomic distribution, we introduce 
the characteristic function of the position (zj) distribu- 
tion defined in terms of its Fourier modes 

1 N 

B n :=-J2 e2mkZ3 - ( 4 ) 

i=i 

Introducing the normalised variable 9, the general defi- 
nition Eq. Q reads, with our normalisation, 

B n = (e 2ine ). (5) 

The direct integration of model , using the same tech- 
niques as outlined in [T^ | (i.e. with thermalization in- 
duced by collisions) , provides the results shown in Fig. ^ 
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Fig. nji shows the appearence of several Fourier modes, 
which reach a steady state amplitude after the backward 
field has grown to its final value (Fig. Qd, solid line). 
We stress that the latter grows from zero spontaneously, 
without any seed being injected in the computation. 

From Fig.^L we remark that one Fourier component of 
the atomic spatial distribution significantly grows away 
from zero well before the backward field begins to de- 
viate from its initial value. We will take this as a sug- 
gestion, to be later verified, that the bunching is respon- 
sible for the appearence of the counterpropagating field. 
The harmonic components, instead, seem to be generated 
by the interaction of the inhomogeneously distributed 
atoms and the optical standing wave generated by the 
two macroscopic counterpropagating fields, (x/ and x b ). 
In addition, we also see how the forward field (dashed 
line in Fig. \\}p) is modified by the presence of the spa- 
tially modulated atomic distribution. 

The physical picture described by this figure is the fol- 
lowing: the cavity initially contains the atoms and no 
fields. At time t = 0, we start injecting the external com- 
ponent Y, collinear with the forward direction. During 
the transient preceding synchronisation, the atoms use 
part of the energy to adapt their positions and velocities, 
thus building up a non homogeneous spatial distribution. 
We remark that the total energy injected into the cav- 
ity (see dot-dashed line) is larger than that contained in 
the two field components at steady state (t > 300 time 
units). The fact that both fields, Xf and x b , grow nearly 
to the same value implies that one cannot neglect the 
dynamics of the forward field. This dynamics maximizes 
the spatial modulation and both fields adapt themselves 
to the global coupling that ensues. It is also crucial to 
remark that, although the spatial atomic inhomogeneity 
appears to be responsible for the generation of the back- 
ward field, at time t w 230 time units, the B\ component 
has reached only one third of its final value. Hence, the 
rest of the modulation (fundamental plus all the harmon- 
ics), results from the interaction with the optical stand- 
ing wave. This point is significant, since it implies that 
the scattering from one field into the other is going to 
be quite symmetric: any increase in one component will 
imply a larger number of photons in one direction, hence 
a larger scattered contribution from that direction into 
the other. The residual difference between the forward 
(dashed line, Fig. QJi) and the backward field is to be 
attributed to the fact that the former receives an addi- 
tional contribution from the external injection; hence the 
corresponding mode contains more photons at any time. 



A more convincing illustration of the active role played 
by the density grating in the generation of the backward 
field, can be gained by decomposing the "force" field act- 
ing on the intensity \x b \ 2 of the backward field into two 
contributions 
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FIG. 1: Fig. (a), shows the time evolution of amplitude of 
the first Fourier modes of the ^-distribution (n = 1,2,3 cor- 
respond to solid, dashed and dotted curves, respectively). 
Fig. (b) shows the rise of both amplitudes (solid and dashed 
line corresponds to \xb\ and \xf\). The dot-dashed line rep- 
resents the injected field, Y. Parameter are: A a = —30, 
T x = 1, r || = 2, C = 7.1, A c = 0, G = 54, Y = 1.45, 
j c = 1/3, T= 5.10" 2 . 



The first term on 



d\x b \ 
dt 



■T{xf,x B ) = -2\x b \ 2 + FA.{xf,Xb), (6) 



where F A (x f ,x b ) = 2C^e((se te )i 
the r.h.s., due to cavity losses, has always a stabilizing 
effect, while F^XfjXb), which accounts for the atomic 
contribution, can be either stabilizing or destabilizing. 
In order to gain some insight on the resulting dynamics, 
we compare the force fields generated in the presence and 
in the absence of bunching, respectively. We do this un- 
der the assumption that the atomic degrees of freedom 
can be neglected j^J, i.e. that the relevant dynamical 
properties are contained in the dependence of Ta on x b 
and Xf. As long as the atomic variables rapidly converge 
to a stationary state, one can treat \xf\ and \x b \ as if 
they were fixed parameters and thereby numerically de- 
termine the "force" field for different choices of the two 
amplitudes. In principle, one should also take care of the 
time-dependence associated with the detuning v between 
the two fields; however, numerical simulations of the full 
model have revealed that the "force" field depends very 
weakly on v. The force resulting for the same parameter 
values as in Figgis plotted in Fig. [21 for different values 
of the backward field amplitude. The solid line with no 
symbols represents the total effective "force" acting on 
the backward field. Steady state operation exists only 
when the resulting total force is zero, and its stability 
is obviously determined by the local slope. It is clear 
that the zero-field state is unstable, while \x b \ w 0.6 is 
a stable one. In this latter regime, one can state that 
the destabilizing action of the atomic contribution, T a 
(dashed curve without symbols in Fig.[2J), is balanced by 
the stabilizing effect due to the cavity loss term, \x b ^ . 

In order to test of the true role played by atomic 
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bunching (i.e., as a master or a slave variable), we have 
proceeded to "switching off' both collisions and optical 
forces acting on the atoms (this is obviously possible only 
in a simulation). The corresponding curves representing 
(dashed line in Fig. [5J and the total force, J 7 , (solid 
line) are identified by symbols. As a result of the absence 
of external perturbations, each atomic position 6j evolves 
linearly according to the velocity value at the "switch off" 
time; in this way, the density grating is rapidly washed 
out. This procedure allows us to separate the possible 
gain component coming from the density grating from 
the one arising in a pumped ensemble of two-level sys- 
tems |22j . The stationary value of J- a under such condi- 
tions remains always negative (so does a fortiori the total 
force) , indicating that the presence of a density grating is 
an indispensable ingredient to maintain a finite backward 
field. Notice that this is at variance with the transition 
studied in [T^], where the force field would be positive 
even in the absence of a density grating [14j. Hence, we 
can safely conclude that the onset of a finite backward 
field is a clear example of CARL and that the steady state 
bunching obtained is the cause for the appearence for the 
backward field component, and not a consequence of the 
existence of an optical standing wave (cf., e.g., 01) m 
the parameter range which we are investigating. 
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FIG. 2: The force Ta (dashed lines) and T (solid lines) versus 
the backward field amplitude x b . Curves with no symbols cor- 
respond to the simulation of the full model; diamonds display 
curves obtained by imposing a flat distribution of positions 6 
(same parameter values as in Fig. . 



IV. LARGE DETUNING LIMIT 

In Section lllll we have just reported purely numerical 
observations. In the perspective of constructing an an- 
alytical description of the transition, it is important to 
identify the possibly few relevant variables. For this pur- 
pose we concentrate on the large detuning limit, where 
the standard adiabatic elimination holds. Although this 
allows us to analytically treat only a limiting case, we will 
draw from it the benefit of later comparing the predic- 
tions of this approach to the recent experimental obser- 
vations [TTj . obtained with very large detuning between 



field and atoms (A a ss 10 6 3> 1). A discussion of the pre- 
dictions resulting from such an analysis will be offered in 
Section El 

By setting, in Eqs. i|T|l. dj = Sj = 0, expanding the 
solutions in powers of 1/ A a (to be considered as a small- 
ness parameter) and considering first order terms, one 
finds dj = — 1 and 

8i = -£(x f e»'+x h e-»>) . (7) 

Accordingly, our model, Eqs. QJ, reduces to the following 
set of equations 

6j = Pj, 

Pi = ^m(x f x* b e 2ie '), (8) 

x f = — [(1 + iA)xf — Y + iCxbB^i] , 

x b = - [(l + iA)x b + iCx f B 1 \ , 

where C = CG/A a , and A = A c + C accounts for the 
combined effect of atomic and cavity frequency shifts in 
the absence of bunching (i.e. B\ = B-\ = 0). The above 
form is similar to the one discussed in |23|. which was, 
however, obtained by phenomenologically adding relax- 
ation terms to a single pass model |24j . 

A. Vlasov Model 

It is now convenient to pass to a Vlasov-like descrip- 
tion by introducing the distribution Q(Q,p) of positions 
and momenta. In order to simplify the analysis we now 
assume that also the atomic positions are randomized in 
each collision |2^|. Accordingly, we can write 

dtQ+pdgQ + ^m(x } xle 2l9 )d p Q = - 7c (Q - Q ), (9) 

where d y denotes the partial derivative with respect to 
y. In the absence of a backward field, Q(0,p) converges 
towards the equilibrium distribution Qq. 

We now study the stability of the solution Xf 
Y ; I + iA) := To, x b = 0, and Q = Qo, by introducing 
the perturbations 6xf, Sx b , SQ. Since the equation for Xf 
depends on the other variables only at second order in the 
perturbations, it can be solved separately. Moreover, dxf 
affects Sx b only at second order, so that it can be factored 
out. Therefore, the stability of Sxf justifies the assump- 
tions made in deriving the single-pass model where Xf is 
considered to be a constant parameter 0, Il2l IT3 | . One is 
thus left with the following two equations 

5x b = -(l + iA)5x b -iCY6B u (10) 
d t 5Q - -d p Q Qm{Ydx* b e 2ie ) - ( lc +pd e )5Q. (11) 

From the structure of such equations we see that if the 
bunching is washed out - i.e., B±i = and SQ = 
no instability can be expected in this re gim e. This is 
at variance with the case discussed in [13, llJI where the 
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effective stability of x b can be turned into an instability 
sufficiently close to the atomic resonance. 

From Eq. one realizes that the coupling occurs 

through the first Fourier mode of the density distribution. 
One can thus solve the above equation by introducting 
the following Ansatze, 



6Q(6,p,t) = Q(p)e 



-2iB+\t 



+ C.C. 



Sx b = £ b e xt , (12 



where A is a complex stability eigenvalue. After substi- 
tuting the above expressions into Eqs. (flOlll II) and equat- 
ing terms with the same exponential factor, one finds 

/>oo 

(A + 1 + iK) £ b + 2mCY / d P Q(p) = 0, (13) 



-Y*£ b d p Q + (X-2ip + 7c )Q = 0. (14) 

By solving Eq. (|14[) and replacing into Eq. (|13|> . one ob- 
tains the solvability condition 



\ = -(X + iA) -irC\Y\ 2 / dp 



dpQo 



A - 2ip + 7 C 



(15) 



Within this framework, the onset of both bunching and 
backward field is signalled by a change of sign of the real 
part of A, A r [2^|. The integral in Eq. (|15fl can be ana- 
lytically expressed using an error function. Nevertheless, 
this would hide an important point. Let place ourselves 
at threshold, and rewrite Eq. (|15|) . assuming that A is 
purely imaginary. Setting A = — we obtain for the 
real part 



dp- 



p e 



-p 2 /(2T) 



tVz^t J-oo (J3 + 2p) 2 + 7c 2 



(16) 



In order for such an equation to be valid, the condition 



dp 



p e 



V/(2T) 



09 + 2pf + 7 2 



> 



(17) 



must be fulfilled. The numerator is an odd function of 
p. The sign of the integral in Eq. (|17|) is thus deter- 
mined by the denominator, whose minimum is reached at 
p = -(3/2. Finally, Eq. JUJ imposes j3 < 0. From the 
sign convention adopted here (the same that in Ref. fifty). 
a complex field E, of frequency to, is linked to its slowly 
varying envelope A as E — Aexp[— iuit]. Thus (3 physi- 
cally represents the detuning between the backward field 
and the forward field. One can thus conclude that the 
spontaneous symmetry breaking causing x b to rise can 
only give rise to a red detuned field. 



B. Fokker-Planck description 

We now turn to analyzing the relaxation mechanism 
that most closely models the recent experimental re- 
sults 11]. There, the cw superposition of additional 



laser beams, tuned to the D2 transition of Rubidium, 
introduces a mechanism which slows down those atoms 
that would be accelerated away from resonance by the 
strong interaction arising from CARL. Indeed, it was re- 
ported [13 that in the absence of the molasses only tran- 
sient backwards emission could take place (and its inter- 
pretation as CARL is complicated by the technique used 
for preparing the sample, which pre-traps the atoms in a 
standing wave). The key point in the experiment is that 
the molasses is always present during the interaction but 
on the other D-line (the experiment being done on the 
Dx). 

In the large detuning limit, the second of Eq. JSJl is 
replaced by 



i>] = %rn(x f x* b e 2l9:i ) - jpj + rjj(t) 



(18) 



Eq. (|18J) describes the evolution of (9, p) in a potential 
self-consistently determined by the dynamics of the two 
counter-propagating field components, Xf and x b , which 
in turn depends on the joint distribution Q(9,p). The 
evolution of Q corresponding to Eq. I|18l) , is given by the 
Fokker-Planck equation 



d t Q + pd e Q + %m{x f x* b e 2t0 )d p Q = 
~fd p { P Q + Td p Q) . 



(19) 



Notice that this equation differs from the Vlasov equation 
(|5J) only for the relaxation terms contained in the right 
hand side. 

By proceeding as in the previous subsection, we lin- 
earize Eq. (|19|) obtaining 



d t 6Q= - d p Qo%m(Y6x* b e 2lf> ) 

-pd e 8Q + jd p ( P SQ + Td p 5Q) , 



(20) 



which is to be compared to the corresponding Eq. 1|11|) 
within the Vlasov description (the linearised x b equation, 
being identical to Eq. ((Tn|) is not reported). Using the 
Ansatze (|12|l . we finally find 



-Y*£ b d p Qo + [X-2ip- 1 )Q 

-l(pd p Q + Td 2 Q) -0 



(21) 



which is to be solved together with the first of Eqs. i|21|l . 
In view of its differential structure, it is doubtful whether 
a general analytic solution can be found, and the devel- 
opment of a numerical analysis would be a quite deli- 
cate task as well. Nevertheless, in the "strong friction" 
(SF) limit, the atomic momentum can be adiabatically 
eliminated and one can thereby perform a quite detailed 
investigation of the transition scenario. While a careful 
discussion of the adiabatic elimination in the presence of 
noise can be found in Ref. [27j , here we limit ourselves to 
presenting a sketchy but substantially correct derivation 
where we set p = in Eq. i|18|) . By proceeding in this 
way, one obtains the Langevin equation 



1 



Xsm{xfX b 



;2,(l 



7 ' 



(22) 
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that is equivalent to the one-variable Fokker- Planck 
(Smoluchowski) equation 



1 



T 



d t p+ -d e %m(x f x* b e 2i9 ) p djp = (23) 



7 

for the variable p 



/oo 
Q(0,P) dp. 
-oo 



(24) 



We remark that Eq. I|22|l . which describes the dynam- 
ics of an ensemble of mean field coupled oscillators, be- 
longs to the class of Kuramoto systems [2^. The main 
differences are: a Dirac-like distribution of eigenfrequen- 
cies and a mean field self-consistently provided by a dy- 
namical equation. By paralleling the approach in Sec- 
tion ^^^J the stability of the solution xj — Y, Xb = 0, 
p = 1/2tt can be studied by introducing the infinitesimal 
perturbations Sxb and S p. The linearisation of Eq. Ij23|l 
yields 



d t Sp 



*„2i0 



27T7 



do^sm (YSxle 



T 

1 



-d 2 e 8p. 



(25) 



Upon then inserting the Ansatze 5p = re 2z0+xt + c.c and 
8xb = £&e A ' into the above equation, one obtains 



£ h Y* 



1 



2tt A 7 + 4T' 



(26) 



and from the linearised equation for the field dynamics, 
Eq. O, 



A = -(1 + iA) + i 



C\Y\ 2 
A 7 + 4T' 



where now, 



(27) 

B±i = I e ±2W p{9) dd. (28) 
Jo 

The threshold can be now obtained by setting A = —i/3 
(where (3 is a real number) in Eqs. I|26I27|) . The resulting 
two real equations for and \Y\ are 

7/3 2 -A 7 /3-4T = 0, (29) 
—(3 (7 + AT) + 4AT = C\Y\ 2 . (30) 

Only the negative (3 solution of Eq. (|29|l leads to a phys- 
ically acceptable solution iFl 2 > 0. The threshold equa- 
tions finally read 



JTH 




(31) 



\Yth\ 



A 1 /A 

^(4T-7) + ^(4T + W - 



T 
4-. 
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The first equation shows that the probe field is always red 
detuned at threshold. This result is in agreement with 



the outcome of the Vlasov model (see Section llV Afl even 
though these two models correspond to different thermal- 
isation mechanisms. 

In the case of a resonant cavity field, A = 0, the above 
equations reduce to 



Pth — — 2^ 



\Yth\ 2 = |(4T + 7)W^. 



(32) 
(33) 



In the SF limit, one can not only determine analytically 
the threshold condition, but can also go beyond, describ- 
ing the regime above threshold. By inserting the expan- 
sion of p, 



(28+pt) 



(34) 



n— — oo 



into Eq. (J21, one obtains the recurrence relation 
7/n = -M/n - n [Rf n _! - R*f n+1 ] - ATn 2 f n , (35) 

where we have set Xfx\ — Re ll3t , R being a complex 
constant. For n > 0, the stationary state is given by 



(4Tn + i 7/ 9) /„ + [Rf n _ x - R*f n+1 ] = . 



(36) 



By introducing Z n = R*(f n+ i/ f n ) and solving the equa- 
tion by means of the continuous fraction method, we ob- 
tain the recurrence relation 



Z n — 1 — 



\R\ : 



Z n - (ATn + i-f/3) ' 



(37) 



The normalisation condition for p gives 27r/ = 1, while 
the atomic contribution for x b and x f can be then deter- 
mined from Eq. 



B 



27r/_xe 



-ipt 



-i/3t 



J,0t Z i0t 

— A* 



(38) 



1 = 2,t/|C " = ~o 



This methods enables us to go beyond the linear analy- 
sis determing the dependence of both Xf and Xb on the 
injected field (see Section[Vj, with the only warning that 
Zq has to be computed numerically. 



V. PREDICTIONS 

In this section we discuss the results derived in Sec- 
tion [W] in the large detuning limit for both the Fokker- 
Planck (see Sec. IIVBII and Vlasov (Sec. IIV A|l models. 
We first test the validity of the two approaches by com- 
paring them to the experimental results of Ref . [Tl] . The 
experiment has been performed with cold 85 Rb atoms 
in a high-Q ring cavity. The atomic parameters are: 
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m = 1.4 10 _29 kg, V = 1.5 lO-^Cm" 1 with relaxation 
rates for atomic polarisation and inversion 2"/± = 711 = 
5.9MHz, respectively. Moreover, the atomic frequency 
is u) = 3.77 x 10 14 Hz (for the D x line), while the in- 
put field is shifted ITHz away from it. The beam di- 
ameter is 130/im. The cavity linewidth is 22kHz, the 
medium's length L = 10~ 3 m, and the transmissivity 
T = 1.8 x 10~ 6 . Since a servo-mechanism continuously 
adapts the dressed cavity resonance to the frequency of 
the injected field, we assume A = 0. 

The equilibrium gas temperature, in the absence of col- 
lective interaction, is experimentally known with a large 
uncertainty. We use an estimated value, f w 250^K, 
which lies within the range of uncertainty. The sec- 
ond crucial parameter which caracterizes the thcrmalisa- 
tion is the damping constant, experimentally evaluated 
in Ref. 01 to be 7 = 9. We will thus use 7 = 7 C — 9. 




0.25 0.5 0.75 1 1.25 1.5 1.75 



P, (W) 

FIG. 3: Bifurcation diagram numerically computed using the 
Fokker-Planck model. Figs, (a), (b), (c), present respectively 
the backward field power, Pb, the forward field power, Pf and 
the detuning between backward and forward field, /3, as func- 
tion of the injected power Pi. All quantities are in physical 
units. The parameters correspond to those of the experiment 

m. 

We first report the results of the Fokker-Planck model 
since it describes more appropriately the molasses dy- 
namics. In Fig. |3K one remarks, from the shape of the 
backward field power, the existence of a continuous phase 
transition. The threshold is located at Pi « 0.2W, com- 
patible with the measurements of Ref. |ll| , which reports 
backward lasing as occurring for a few Watt of injected 
power. 

The possible underestimate in the threshold value, 
which we obtain from our approach, can be attributed 
to two experimental features that cannot be taken into 
account in a simple way. First, the servo-control which 
keeps the cavity on resonance by maximizing the forward 
transmitted field [Tl| cannot be correctly modeled by sim- 
ply setting A = 0. Indeed, its operation amounts to 
a dynamical process which counteracts the birth of the 
backward emission, since the latter removes energy from 
the forward field. As a result, one expects the thresh- 



old to be displaced towards higher input power values. 
Second, the spatial modulation imprinted onto the sam- 
ple by the molasses beam, with a wavelength that differs 
from that of the collective process (by a factor close to 
V2) , introduces a competition between two incompatible 
length scales. It is likely that this effect should also con- 
tribute to raising the critical input power even further. 

The detuning value, /3, at which the backward field 
arises (cf. Fig. also agrees with the observations 0] , 
which report a red-shift around 100kHz. Moreover, the 
continued fraction expansion discussed in the previous 
section also allows for a determination of the behavior 
of the forward field as a function of the input power (cf. 
Fig-Eb)- Below threshold, Pp grows linearly as a func- 
tion of the injected power, Pj, as no energy is transferred 
to the backward field. Above threshold, a fraction of the 
injected field feeds the backward mode, thus leading to a 
slower growth of Pp. 




Pj (W) 

FIG. 4: Bifurcation diagram numerically computed using the 
Vlasov model. Figs, (a), (b), (c), present respectively the 
backward field power, Pb, the forward field power, Pf, and 
the detuning between backward and forward field, /3, as func- 
tion of the injected power Pi. All quantities are in physical 
units. The parameters correspond to those of the experiment 

m. 

Although one does not expect the Vlasov description 
to closely reproduce the experimental results reported in 
Ref. |ll|. it is worth analyzing its predictions; this will 
help us clarifying the role of the relaxation mechanisms 
and will sharpen our understanding of the problem. From 
Fig. 0] one sees that, for the same parameter values used 
in the Fokker-Planck approach, the threshold is some- 
what lower and the red shift in the backward field fre- 
quency is smaller as well. Nevertheless, the differences 
are not so dramatic, since, while the threshold occurs for 
somewhat lower values of the input power, the frequency 
shift shows an even somewhat better agreement with the 
observations [TT| . In addition, the uncertainty on some 
of the experimental parameters is large enough to render 
a sure discrimination between the modeling based on a 
Vlasov or on a Fokker-Planck process quite difficult. Be- 
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cause of the intrinsic paradox, given by the inclusion in 
the experiment of the optical molasses which ought to 
exclude the Vlasov description, this result strongly indi- 
cates that the current state of the system's characteriza- 
tion, both experimental and theoretical, is not adequate; 
hence, a large amount of caution in drawing conclusions 
ought to be exercised. 



0.2 
0.1 



P=0.25VK 





) 


2 


4 6 






l 

„ Pj=1.5W_ 












i , f 



0.2 
0.1 





' 1 1 1 
/\ P=2W /\_ 








, 1 I f 




e 



FIG. 5: Position distribution p(6), for different values of the 
injected power Pi, in physical units. In each part of the fig- 
ure, the dotted line corresponds to the uniform distribution. 
Parameters correspond to the experimental values of Ref. [Til ]. 



above threshold is plotted in Fig. |BJ where one can in- 
deed appreciate the need to account for both the p and 
9 dependence in this framework. Nevertheless, one can 
compute p(9) by integrating Q(9,p) (see Eq. {23J0- The 
results, plotted in Fig.0 reveal again a qualitative agree- 
ment with the Fokker-Planck approach. The larger mod- 
ulation for the same input power is a consequence of the 
smaller value of the threshold in the Vlasov description. 




FIG. 7: Atomic position distribution p(9), as function of 6, 
above threshold. The dotted line corresponds to the uniform 
distribution, below threshold. 



In addition to the direct comparison to the experiment, 
we can also examine some specific predictions coming 
from each model. The distribution of atomic positions p 
predicted by the Fokker-Planck model is plotted in Fig.|5] 
the different parts of the figure correspond to increasing 
values of the injected power Pj (the other parameters 
remaining unchanged). One can clearly see the rise of 
the modulation with an increasing role played by higher 
harmonics. 

In the Vlasov description, there is no equivalent regime 
characterized exclusively by the distribution of 9- values. 
The probability density Q(9,p) for an input intensity 




FIG. 6: Joint distribution Q(9,p), above threshold. The 
rameter values are the same as those of the experiment | 
with input power of 480mW. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we have presented a general framework 
for the description of the spontaneous appearence of a 
stationary density grating arising from the interaction be- 
tween an atomic sample enclosed in a bidirectional ring 
cavity and an injected, quasi-resonant field. This con- 
tribution completes the picture previously sketched out 
m Refs. [U H3, where we have shown that the inclu- 
sion of a relaxation mechanism for the external degrees 
of freedom is indispensable for a correct modeling of the 
interaction, and is responsible for the appearence of a 
stationary backward field. In the previous work, we fo- 
cused on the existence of a different kind of phase tran- 
sition, initiated by a grating in the phase of the atomic 
polarization. Here, with the help of a detailed numerical 
analysis we have proven beyond doubt, that a true sta- 
tionary CARL action - i.e., initiated by the birth of an 
atomic density modulation - can spontaneously occur. 
For the investigation, we have chosen to study a bidirec- 
tional ring cavity model 0, because of the fuller and 
more correct description that the latter provides of a real 
system. In the course of the presentation, we have shown 
that in certain limits the ring cavity model reduces to the 
single-pass one 0,0,0; this legitimizes, a posteriori, 
the use of the latter approach under such conditions. 

Following a discussion of the main relaxation mecha- 
nisms that one may expect to find in an experiment, we 
have offered an analytical description of the two most 
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likely ones on the basis of a probabilistic description, 
in the large detuning limit. Comparison between our 
predictions and recent experimental results 0] shows a 
semi-quantitative agreement, surprisingly, for both the 
Fokker-Planck and the Vlasov approach. Besides the 
need for a more accurate experimental characterization 
of the phenomenon, including a better determination of 
the parameter values, and systematic measurements of 
the threshold behavior and of the functional dependence 
of the backward field power and frequency, additional 
theoretical work needs to be completed. In particular, a 
more accurate modeling of the relaxation mechanism rep- 
resented by the optical molasses needs to be achieved, to- 
gether with an analysis of the limits of validity of the de- 
scriptions proposed for the losses. Finally, it will be use- 
ful to investigate the transtion scenario when two mech- 



anisms separately described here and in Refs. 
cooperate to induce a backward field; this is likely to oc- 
cur in the range of smaller detunings when the atomic 
variables cannot be eliminated. 
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